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ON  THE  DYNAMICS  OF  FLOATING  FOUR-BAR  LINKAGES  * 


Rui  Yang 


P.  S.  Krishnaprasad 


Electrical  Engineering  Department 

& 


Systems  Research  Center 
University  of  Maryland,  College  Park. 


ABSTRACT 

The  hamiltonian  structure  of  floating ,  planar  four-bar  linkages  is 
discussed.  The  geometry  of  configuration  space  is  related  to  the  classical 
theory  of  mechanisms  due  to  Grashof.  For  generic  value  of  kinematic 
parameters,  the  techniques  of  symplectic  (and  Poisson)  reduction  apply. 


1.  INTRODUCTION 

There  has  been  significant  progress  in  our  understanding  of  the  hamiltonian  structure 
of  serial-link  (or  open  chain)  multibody  systems  [2,3,9,11,12,15-18,21-23].  The  use  of 
geometric  methods,  symmetry  principles  and  reduction  has  led  to  deeper  knowledge  of  the 
phase  portraits  of  model  problems.  This  insight  has  been  helpful  in  developing  appropriate 
control-theoretic  tools.  The  primary  source  of  motivation  for  these  problems  has  been 

*  This  work  was  supported  in  part  by  the  AFOSR  University  Research  Initiative  Pro¬ 
gram  under  grant  AFOSR-87-0073  and  by  the  National  Science  foundation’s  Engineering 
Research  Center  Program:  NSFD  CDR  8803012. 
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in  aerospace  engineering  where  imaginative  designs  of  multibody  spacecraft  have  been 
proposed  and  on  occasion  realized  [4,5,24]. 

In  contrast,  multibody  systems  with  kinematic  loops  present  serious  challenges.  The 
loop  constraint  may  lead  to  singular  configuration  spaces.  The  Dirac  theory  of  constraints 
[6,20]  applies  in  the  smooth  setting.  Little  or  nothing  is  known  about  hamiltonian  structure 
and  phase  portraits  in  concrete  cases.  However,  engineering  applications  suggest  that 
multibody  systems  with  kinematic  loops  are  of  practical  importance  [10].  Parallel  linkage 
based  robot  manipulators  are  contemplated  for  space  applications. 

Here  we  discuss  the  geometry  and  dynamics  of  floating  four-bar  linkages.  In  this 
model  problem,  the  classical  Grashof  criterion  [8,10,14]  appears  through  conditions  for  the 
configuration  space  to  be  a  smooth  manifold.  The  topology  of  the  configuration  space  is 
also  related  to  the  Grashof  criteria.  We  explore  symmetry  properties,  hamiltonian  structure 
and  reduction  of  four  bar  linkage  dynamics.  Explicit  computation  of  constrained  dynamics 
is -difficult.  Yet  in  the  present  setting,  using  geometric  techniques,  one  can  infer  qualitative 
properties  without  recourse  to  explicit  analytic  representation  of  the  constrained  dynamics. 
We  use  a  theorem  of  Smale  to  determine  relative  equilibria  for  four-bar  linkage  dynamics. 


2.  NOTATIONS  fe  GEOMETRIC  CONSTRAINTS 

The  structure  of  a  closed  floating  four-bar  linkage  is  represented  in  Fig.  1.  The  bars 

are  labeled  clockwise  from  0  to  3  as  shown.  We  define  the  following  quantities. 

dty  the  vector  of  hinge  point  which  connects  i-th  bar  with  j-th  bar 

relative  to  a  body-fixed  frame  with  origin  at  the  center  of  mass  of 
the  i-th  body; 

r,  the  position  vector  of  the  center  of  mass  of  i-th  bar  relative  to  an 

inertial  observer; 

rf  the  vector  from  the  system  center  of  mass  to  the  center  of  mass  of 

z  -th  bar; 
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inertial 

observer 


Fig.  1  The  general  structure  of  four-bar  linkage 


rc 


R(9i) 


the  position  of  the  system  center  of  mass  relative  to  the  reference 
point  of  the  inertial  observer; 


the  rotation  through  angle  9i  giving  the  orientation  of  i-th  bar 


relative  to  the  inertial  space; 


R(0i)  = 


(cos(9i) 
sin(6i ) 


cos(6i)  )  ' 


R($ij)  joint  rotation  between  i-th  and  j-th  bar, 


R(9tj)  =  R(9i  -  9j)  =  R{9i)R{-9j)- 


u 


the  distance  between  the  joints  on  i-th  bar,  (or  ’’length”  of  i-th  bar); 
all  5”  0 ;  1%  —  Hd^i+i  , 
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Ii  the  moment  of  inertia  of  i-th  bar  about  its  center  of  mass; 

m  the  total  mass  of  the  system,  i.e. 

3 

m  =  m,-. 

2  —  0 

With  the  above  notations,  any  pair  of  adjacent  bodies  is  connected  by  the  following 
relation,  so  called  hinge  constraint, 

r-+1  =  ri  +  iZ(0i)di)i+i  -  i  =  0, 1,2,3  (mod  4).  (2.1) 

By  eliminating  r?  in  Eq.(2.1)  we  find  the  loop  constraint  or  closure  constraint , 

3 

Y  W)(dM-+i  -  dM-i)  =  0,  (2.2) 

t=0 

where  we  adopt  the  convention  that  d,4  =  d,-,o  and  d4>,-  =  do,,- . 

3.  THE  CONFIGURATION  SPACE 

In  this  section  we  investigate  the  conditions  under  which  the  loop  constraint,  Eq.(2.2), 
describes  a  submanifold  with  respect  to  the  configuration  manifold  of  an  open  four-bar 
chain. 

For  a  planar  floating  four-bar  open  chain,  the  configuration  space  is 

M  =  R2  x  S1  x  S1  x  S1  x  S\ 

where  M  is  a  6-dimensional  smooth  manifold  with  local  coordinates  of  the  form, 

q  =  (so, i/o,0o, 0i, 02, 03). 

This  corresponds  to  keeping  track  of  a  material  point  (say  center  of  mass  on  one  of  the 
bodies)  and  the  four  absolute  orientations.  See  [12,21,22]  for  the  hamiltonian  mechanics 
of  such  open  chains. 
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For  a  closed  four  bar  mechanism  as  considered  in  this  paper,  the  configuration  space 
denoted  by  Q  is  a  subset  of  M  determined  by  Eq.(2.2),  or  simply 

3 

Q  =  {q  €  M|F(q)  =  mx*i,i+ 1  -  difi- 1)  =  0}.  (3.1) 

t=0 

Note  that  F  :  M  — >  R 2  ,  and  from  [i]  we  know  that  if  0  is  a  regular  value  of  F,  i.e. 
dF/dq  has  full  rank  for  all  q  6  Q,  then  Q  is  a  submanifold  of  M . 

From  (3.1)  we  have, 

_  ( 0  0  —l0sin(6o)  — lisin(8i )  — l2sin(02 )  —hsin(03)\ 

dq  V°  0  locos(80)  licos(8i )  l2cos(d2)  hcos(83)  J  ' 

Then  it  is  easy  to  check  that  all  the  nontrivial  determinants  of  2  x  2  submatrices  are  given 
by  the  set  of  expressions, 

fl'i(q)  =  khsin(8i  -  80) 

02 (q)  =  lohsin(82  -  80) 

03(q)  =  lohsin(83  -  80 )  (3.2) 

04(q)  =  hhainfa  -  ^i) 

05 (q)  =  hhsin($z  -  0i) 

06 (q)  =  hhsin(03  -  Q2). 

Therefore  for  each  q  G  Q ,  if  there  exists  an  i  such  that  gi(q)  ^  0  ,  Q  is  a  submanifold  of 
M. 

It  is  obvious  that  the  above  condition  depends  on  the  lengths  of  the  links.  To  find 
the  condition  on  the  links,  we  can  use  an  equivalent  way,  that  is,  find  necessary  conditions 
on  the  lengths  of  the  links  such  that  <7,-(q)  =  0  for  all  i . 

From  Eqs.(3.2)  it  can  be  seen  that  if 

8i  ~  Oo  =  0  or  7r  (3.3)a 

and 

83  —  82  =  0  or  tt  (3.3)6 
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and 

03  -  #0  =  o  or  7T  (3.3)c 

then  <?i(q)  =  0  for  all  i. 

Premultiplying  Eq.(2.2)  by  7?,((90) ,  we  get  following  equivalent  closure  constraint 
equations 

Iq  +  l\cos(6i  -  90)  +  h(cos(9i  -  92)cos(03  -  0o)  +  sin(03  -  82)sin(63  -  60)) 

+l3cos(63  -  62)  =  0  (3.4 )a 

lisin(8i  —  #0)  +  h  (cos(93  —  82)81x1(83  —  90)  —  sin{93  —  82)cos{83  —  8q)) 

+l3sin(83  -  82)  =  0  (3.4)6 

The  conditions  given  in  Eq.(3.3)  make  Eq.(3.4)6  trivial.  By  choosing  all  the  possible  value 
of  8i—80,  83  —  82  and  83  —  80  given  in  Eq.(3.3),  from  Eq.(3.4)a,  we  can  find  the  link 
conditions,  which  are  summarized  in  Table  1. 

Theorem  3.1:  If  /0  dh  l\  ±  h  ±  h  7^  0,  Q  is  a  submanifold  of  M . 

From  Table  1,  it  is  easy  to  observe  that  case  (i)  can  never  happen  since  li  are  assumed 
to  be  positive.  In  addition,  cases  (iii),  (iv),  (v)  and  (vi)  are  trivial  cases  since  none  of  them 
can  be  formed  by  any  general  four-bar  closed  loop.  In  these  cases  the  configuration  spaces 
lose  one  degree  of  freedom  and  are  three  dimensional. 

Furthermore,  Theorem  3.1  can  be  simplified  by  ignoring  the  labels  on  the  bars.  To 
do  this,  we  first  recall  some  definitions  and  results  in  the  classical  theory  of  mechanisms 
[8,13,14].  We  define  following  quantities 

s  =  length  of  shortest  bar 
l  =  length  of  longest  bar 
p,q  =  lengths  of  intermediate  bars. 

A  link, which  is  free  to  rotate  through  27r  with  respect  to  a  second  link  is  said  to 
revolve  relative  to  the  second  bar  and  is  referred  to  as  a  crank.  Any  bar  which  does  not 
revolve  is  called  a  rocker.  If  it  is  possible  for  all  bars  to  become  simultaneously  aligned, 
such  a  state  is  called  a  change  point  and  the  linkage  is  said  to  be  a  change-point  mechanism. 
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Table  1. 


case 

0}—  0q  0-j—  ©2  0y  0q 

link  condition 

i 

0  0  0 

1  0+l  l+l2+ 13~° 

i  i 

0  0k 

lO  +  ll~l2'l2=0 

i  i  i 

0  k  0 

lO+ll'l2+l3=0 

i  V 

0  K  K 

l  Q+  l  J  +  1 2  -  /  0 

V 

K  0  0 

l  o"  l  Z  2+  ^  3=  0 

v  i 

K  0  K 

l0'll'l2l  f° 

vii 

n  n  0 

Iq-1  pl  tflfzO 

viii 

7C  K  71 

irh+irifO 

structure 


0  n  1 

Q  „  n 

n 5  2 

3  n  0  n  1 
CEZ  ~  . n 


On  1  n  2 

tr . vw  . u - n 


2  n  3  n  0 

a . . . .  n 


1  n  2  n  3 

(T - . n 

0 


3  n  0 

ct — . 1) 


n 3 - r 

»  D 


^  t }  2 _ — 


or 


0 


j±bdrp 


D 


Theorem  3.2 i(Grashof)  (1)  A  four-bar  mechanism  has  at  least  one  revolving  link  if 

s  +  /  <  p  +  q 


and  all  three  will  rock  if 


s  -f  l  >  p+  q- 

(2)  A  four-bar  mechanism  is  a  change-point  mechanism  iff 

s  +  /  =  p  +  q. 

Remark:  It  is  easy  to  check  that  the  cases  (ii),  (vii)  and  (viii)  in  Table  1.  correspond  to 
s  +  l  =  p  +  q  i.e.  they  coi'respond  to  change-point  mechanisms. 
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Corollary  3.3:  If  /  <  s  +  q  +  p  and  s  +  l  ^  p  +  q,  Q  is  a  submanifold  of  M . 

Remark:  Note  that  if  l  >  s  +  q  +  p,  the  mechanism  is  not  constructible. 

In  order  to  find  a  topological  description  of  the  configuration  manifold  Q,  we  first 
introduce  following  fact. 

Proposition  3.4:  If  s  +  /  <  p  -f  q  and  l\  =  s,  i.e.  h  —  Taino<i<3(h) ,  then  93  —  92  kn 
for  k  €  Z . 

Proof.  The  mechanism  can  be  assembled  with  s  adjacent  to  l  or  with  s  opposite  l.  And, 
l  can  be  /o,  h  or  I3  .  If  63  —9?  =  kn ,  the  whole  structure  attains  a  triangular  shape  which 
has  the  property  that  the  sum  of  two  sides  is  larger  than  the  third  one.  Then  it  is  easy  to 
check  that  all  possible  cases  will  lead  to 

s  +  l  >  p  +  q. 


This  is  a  contradiction. 

i 

From  Grashof’s  theorem  and  the  above  proposition  we  can  get  a  topological  descrip¬ 
tion  for  Q . 

Corollary  3.5:  For  a  crank  mechanism  with  s  +  l  <  p  +  q,  each  connected  component  is 
diffeomorphic  to 

Q  =  R2  x  S1  x  S1 

with  parametrization  ( xo,yo,9o,9i )  if  /1  =  s  is  assumed. 

Proof.  We  just  need  to  prove  that  9 2  and  #3  can  be  uniquely  determined  by  90  and  . 
The  constraints  given  in  Eq.(3.1)  give  the  relations  among  9{ ’s  and  they  have  continuous 
partial  derivatives  with  respect  to  #,■ .  Since 

OF 

det(d(02  ))  =  l2hsin(03  -  92), 

it  follows  from  the  Implicit  Function  Theorem  that  if 

$3  —  92  7^  ki r  for  k  G  Z 

in  some  neighborhood  of  a  point  (00, 9i  ,92,93),  then  there  exists  a  unique  pair  of  functions 
/1  and  f2  such  that 

92  =  fi(0o,9i)  and  93  = 


S 


and 

F(£MiJi(0o,0iU2(0o,0i))  =  O. 

Using  the  result  in  above  proposition,  the  proof  is  complete.  In  addition,  as  seen  from  [13], 
there  exists  a  pair  of  function  fi  and  such  that 

#2  =  +  /i(0i  —  #o)  and  O3  —  6q  +  /2(0i  —  Oq).  (3-5) 

I 

Remark:  Following  Gibson  and  Newstead  [7]  one  can  show  that  if  l  +  s  <  p  +  q,  then 
there  are  two  connected  components  of  Q  and  if  s  +  l  >  p  +  q,  the  space  Q  is  connected 
and  diffeomorphic  to  R2  x  S'1  x  S’1 . 

Since  the  configuration  space  under  the  condition  given  in  Corollary  3.5  has  an 
explicit  parametrization,  in  the  rest  of  this  paper  we  shall  study  the  system  on  this  space 
mainly. 


4.  KINETIC  ENERGY 


In  this  section  we  shall  derive  the  kinetic  energy,  or  Lagrangian  since  we  assumed 
that  no  potential  energy  is  involved,  for  the  whole  system.  The  basic  idea  is  to  write 
the  kinetic  energy  for  each  individual  body  first  and  then  use  the  constraint  equations  to 
eliminate  extra  variables. 

The  kinetic  energy  of  the  i  -th  bar  is 

Ti  =  \  +  imillr.-H2 


where  =  Oi .  The  total  kinetic  energy  is 

3 

I  1 

T 

V.  *- 

=0 


i=0  i=0 


(4.1) 


To  describe  the  kinetic  energy  relative  to  the  center  of  mass,  we  have  following  useful 
equations, 


r,  =  rc  +  r  \ 
3 

m i =  0. 

!  =  0 
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i  =  0,1, 2, 3 


(4.2) 

(4.3) 


(4.4) 


By  applying  Eqs.(4.2)  and  (4.3),  Eq.(4.1)  becomes, 


Applying  Eq.(2.1)  and  Eq.(4.2),  we  get 

r-  =  — [.R(0i_i)mi_.1d;_i  * 
m 

-  i?(^*)(mi_i d,,*—!  +  (mi+i  +  mi+2)d{>l+1) 

)((7Tli^-i  *j"  ™-i+2)&i+\,i  ^i+2^i-\-l  ,i-{-2^ 

+  -R(#;+2)?™i+2d»+2,j+l] 

for  i  —  0, 1,2,3  (mod  4).  Furthermore, 

—  d?,/2(0j)(mj_i  dqj-i  +  (m,_|_i  +  m,,+2)di,i+i) 

+  u>i+ii2(^a+i)((mJ+i  +  m,+2)di+iit-  —  m,-+2d,-+i,i+2) 

+  uii+2R(0i+2)mi+2di+2,i+i]  (4.5) 

for  i  =  0,1, 2, 3  (mod  4),  where 


By  substituting  the  formula  for  rf  into  Eq.(4.4),  we  get  a  more  compact  expression  for 
the  kinetic  energy 


T  =  i  <  w,  Jd>  >  +^m||rc||2  (4.6) 

where  Co  =  (wo,o;i,a;2,u>3)T  and  J  =  ( =  0, 1,2,3)  is  a  4  x  4  symmetric  matrix, 
with  elements  given  as  follows. 

Let 


M/J  = 
M(n  = 


—  [m,m!+i(mj  -f  mf+1) 
mz 

+  +  rriimi+ 2) 

+  mim,+2(mt+i  +  mi+2)] 

2  x 

^-(mi+2  -  m.+xm,-!) 
mjm,+2 


ml 


L(mi+i  +  m,-i). 
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(4.7) o 

(4.7) 6 

(4.7) o 


Then 


Ju  =  Ii  +  M{  ||d||?  t+1  +  M/_j  ||d||?  t-_j 

-  2 M/4  <  >  (4.8)a 

J%, i+l  —  -^t  ^  > 

+  M-1  <  d,- ^+1,  -R(0H-i,i)d,- +i,i+2  > 

+  <  diii_i,i?(^i+lj,)d,'4-ijj  > 

+  M{n  <  dj(;_i, i?(0i+x,i)di+i)j+2  >  (4.8)6 

Ji,i+2  =  —Ml1  <  di)J+i,i?(^+2,i)di+2,i+l  > 

Mj_ (.j  <  dj)t-4-i ,  R(6i+2,i)di+2,i-l  ^ 

—  Mil 2  <  di,«-l>-R(^»+2,i)d*'+2)i-l  > 

—  Ml11  <  diti-i,R(6i+2,i)di+2,i+i  >  (4.8)c 


for  i  =  0, 1, 2, 3  (mod  4). 

Up  to  now,  we  have  not  applied  the  condition  on  the  lengths  of  links  of  a  closed  loop. 
It  is  clear  that,  in  Eq.(4.6),  8i  and  the  velocities  u>j,  i  =  0,1, 2, 3  are  involved.  Since  the 
chain  is  closed,  these  variables  are  not  independent.  As  shown  in  the  proof  of  Corollary  3.5, 
if  s  +  l  <  p  +  q  and  l\  =  s,  we  have  an  one-to-one  map  from  (80,8i)  to  (#2, 63).  Therefore 
every  element  of  the  matrix  J  can  be  expressed  as  a  function  of  81  —  60  uniquely. 

Under  the  same  conditions  on  links,  i.e.  s  +  l  <  p  +  q,  the  loop  constraints  (2.2) 
yield  a  relation  between  (wo,wi)  and  (u>2,u>3): 

(:)"<:) 

where 


(l0sin(6s—8o)  _hsin(03—6i)  \ 
l^sinl$3—62)  l2sin(63—02)  | 

losin($2—0 0)  l\sin{82  —  #i)  I  ' 

/3sm(03— $2)  l3sin{83  —  $2)  / 


(4.9)6 


Here,  the  matrix  0  is  well  defined  because  of  the  result  in  Proposition  3.4.  Again,  since 


s  +  l  <  p  +  q  and  s  =  1 1 ,  elements  of  matrix  Q,  are  functions  of  81  —  80 . 


We  summarize  the  above  discussion  in  the  following  theorem. 
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Theorem  4.1:  If  s  +  /  <  p  +  q  and  l\  =  5,  the  kinetic  energy,  or  Lagrangian,  can  be 
represented  as 

T  =  L  =  \  <  >  +^ml|rc|j2  (4-10) 

where  w  =  (wo^i)71  and 

j  =  (/  sF)j(') 

for  J  given  in  (4.6)-(4.8)  and  ft  given  in  (4.9).  In  addition,  the  elements  of  J  are  the 
functions  of  (9i  —  9o). 

Before  ending  this  section,  we  give  a  property  of  the  matrix  ft  which  will  be  used  in 
section  6. 

Proposition  4.2:  Under  the  assumptions  of  Theorem  4.1, 


Proof:  Premultiplying  Eq.(2.2)  by  R(6 2)  and  R{6 3),  we  get 


losin(0Q  —  9 2 )  +  hsin($i  —  02 )  +  l2sin{92  —  92)  —  0 


and 

l0sin(90  -  92)  +  l1sin(9i  -  93 )  +  l2sin(92  -  03)  =  0 
respectively.  They  immediately  imply  the  claim. 
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5.  SYMMETRY  AND  INTEGRAL 

We  shall  show  here  that  a  floating  four-bar  linkage  is  a  simple  mechanical  system 
with  symmetry  in  the  sense  of  Smale  [1,19]. 

A  simple  mechanical  system  with  symmetry  is  a  4-tuple  ( Q ,  K,  V ,  G ) ,  where, 

(i)  (Q,K)  is  a  Riemannian  configuration  manifold  with  metric  K ; 

(ii)  G  is  a  Lie  group  acting  on  Q  on  the  left, 

(g,q)  h->  $g(q)  =  $(g,q) 
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such  that  for  each  g  6  G ,  is  an  isometry  of  ( Q,K ); 

(iii)  V  :  Q  — *  R  is  a  G-invariant  potential  function. 

The  associated  Lagrangian  is  defined  by 

L  :  TQ  — ►  R 

1  (5.1) 

Vg  L(vq)  =  -Jf(ug,v9)  -Vo  r(v9) 

where  r  :  TQ  -»  Q  is  the  canonical  tangent  projection.  The  Legendre  transform  FL  of  L 
is  given  here  by  the  vector  bundle  isomorphism 

Kb  :  TQ  -»  T*Q 


satisfying 

Kb(vq)  ■  Wg  =  K(vq,Wq)  VVg,  Wg  E  TgQ. 

The  hamiltonian  H  :  T*Q  ->  R  is  defined  to  be 

#(«,)=  ljtf((jrk)-»(o,),(Kk)-1(«f))  +  VoT'(o,)  (5.2) 

where  r*  :  T*Q  — >  Q  is  the  canonical  cotangent  projection.  Letting  denote  the 
canonical  symplectic  structure  [l]  on  T*Q,  the  hamiltonian  dynamics  on  T*Q  is  given  by 
the  unique  vector  field  Xfj  on  T*Q  such  that 

dH(Y)  =  tl0(XH,Y) 

for  all  vector  field  Y  on  T*Q . 

Let  :  G  x  T*Q  — *  T*Q  be  the  cotangent  lift  of  the  action  $.  Denoting 

= we  have 

H  =  H,  (5.3) 

i.e.  the  group  G  is  a  symmetry  group  of  the  hamiltonian  system  (T*Q,Qo,Xh)- 

The  modern  setting  of  Nother’s  theorem  relating  symmetry  to  the  existence  of 
integrals  of  motion  is  given  by  the  concept  of  momentum  mapping.  Let  Q  denote  the 
Lie  algebra  of  G  and  Q*  its  dual.  The  map 

J  :  T*Q  — +  Q* 
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given  by 


J(«,)  ■«  =  «,•  ?«(s)  v«,  €  r;a 


(5.4) 


is  Ad*  -equivariant,  where  £  E  G  and  £q  is  the  infinitesimal  generator  of  $  on  Q  associated 
to  £  (see  [x]  corollary  4.2.11).  J  is  a  momentum  mapping  and  the  G-invariance  of  the 
hamiltonian  H  (Eq.(5.3))  implies  that  J  is  an  integral,  i.e.  it  is  conserved  along  trajectories 
of  XH- 

Returning  to  four-bar  linkages,  we  restrict  attention  for  the  moment  to  linkages  of 
Grashof  type,  i.e.  the  condition  s  +  /  <  p  -f  q  holds.  To  fix  the  parameterization,  we  let 
li  =  s.  We  also  restrict  ourselves  to  a  connected  component  of  the  configuration  space. 
Then 


Q  =  R2  x  S1  x.  S1; 


the  kinetic  energy  metric  is  K  defined  by 

K((v,w),(v,w ))  =  im  <  v,v  >  <  w,Jw  > 

(see  Eq.(4.10)  of  the  previous  section).  The  group  of  rigid  motions  in  the  plane  is  the 
symmetry  group: 

G  =  S1  x  R2. 


Denoting  a  point  in  Q  by  (rc,0o,#i),  the  action  $  of  G  =  S1  x  R2  on  Q  is  given 
by 

$((^,r),(rc,0o,0i))  =  (r  +  rc,0o  +  <Mi  +  <f>)- 

For  our  purposes  it  is  convenient  to  eliminate  the  effect  of  translations  altogether  by  putting 
the  inertial  observer  at  the  center  of  mass  of  the  system,  i.e.  rc  =  0.  In  [22]  this  process  is 
explained  via  symplectic  reduction  by  the  translation  group  R? .  The  effect  of  taking  this 
step  is  that  now, 


and, 


0  =  S1x51,  G  =  S1  (5.5) 

$(<£,(0o,0i))  =  (0o  +  <Mi  (5-6) 
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The  kinetic  energy  metric  on  Q  is 


K(wi,w2 )  =<  w1,Jw2  > 


(5.7) 


for  w i,  w2  €  T9Q. 

The  hamiltonian  is  given  by 

H  =  -<(*,  J“V  >  (5.8) 

for  ^  =  K^(w).  Using  the  abstract  formula  (5.4)  and  the  action  (5.6),  one  can  show  that 
a  momentum  mapping  for  the  action  $  is 

J(9,n)  =  no  +  A*i-  (5.9) 

Of  course,  v  =  fi0  +  Vi  is  conserved  along  trajectories  of  Xh  for  the  hamiltonian  H  in 
(5.8)  and  it  is  simply  the  net  angular  momentum  of  the  floating  four  bar  linkage  relative 
to  an  observer  at  the  system  center  of  mass. 

The  dynamical  trajectories  are  confined  to  level  set  of  the  form  The  group 

S1  viewed  as  the  isotropy  subgroup  of  the  momentum  value  v ,  acts  freely  on 
and  one  gets  the  symplectically  reduces  dynamics  Xjjv  on  the  reduced  phase  space 
P„  =  J~l{u)/Sl  ~  51  x  R1.  We  discuss  this  further  in  the  next  section. 


6.  f REDUCED!  DYNAMICS  fe  RELATIVE  EQUILIBRIA 

As  in  [22]  it  is  possible  to  Poisson-reduce  the  dynamics.  We  recall  that  given  a 
symplectic  manifold  (M,w),  and  a  smooth,  free,  proper,  symplectic  action  of  Lie  group  G 
on  M ,  the  canonical  Poisson  structure  on  M  defined  by 

{Mm  =  u <(Xf,X,)  V/,  g  €  C°°(M) 

descends  to  a  Poisson  structure  on  the  quotient  P  =  M/G.  The  latter  is  defined  by 

{f,g}M/G0v  =  {f0n,g°n}M  (6.1) 
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where/, g  €  C°°(M/G )  and  tc  :  M  — ►  M/G  is  the  canonical  projection.  If  H  :  M  — ►  R  is 
a  G'-invariant  hamiltonian,  it  induces  H  :  M/G  — ■*  R  defined  by  H  o  tt  =  H .  We  refer  to 
the  dynamics  (vector  field)  X n  defined  by 


Xfi{f)  =  {f,H}  V/  €  C°°(M/G)  (6.2) 


as  the  Poisson  reduction  of  the  dynamics  Xh  . 

In  the  present  context,  with  Q,K,H  as  in  (5.5)-(5.8),  the  space  M  —  T*(Sl  x  S’1) 
with  parameterization  (9q,8i,  /j,q,  hi)  carries  the  Poisson  structure, 


if  >9)  = 


?— 0 


df  dg  df  dg 
ddi  dm  dm  ddi 


) 


(6.3) 


for  all  /,  g  G  G°°(T*(S1  x  S1)).  The  action  of  G  =  S1  on  Q  given  by  (5.6)  is  free  and 
proper.  The  quotient  P  =  T*(S1  x  S1)/S1  ~  S1  x  R2  carries  a  reduced  Poisson  structure. 
Parameterizing  P  =  T*(S1  x  S1)/Sl  by  8XQ  =  (#i  -  80,/j, 0,  A*i) ,  the  Poisson  bracket  on  P 
is  given  by, 


{/>  9} 


T'^S'xS^/S1 


df  dg  dg  dg  df  df 

d6io  [dm  dfi0  d910  ^ dm 


(6.4) 


which  is  a  noncanonical  structure.  The  reduced  hamiltonian  H  is  given  by 


£(Wo,/*i)  =  H(poA,nhl*i), 


(6.5) 


since  the  matrix  J  in  (5.8)  is  a  function  of  the  difference  #10  =  —  &o  only.  The  reduced 

dynamics  is  then  immediately  given: 


m  - 


dH 
98 10 

dH 


8 


98 10 

dH  dH 


(6.6) 


10 


dm  9 no ' 

Equation  (6.6)  involves  complicated  analytic  expression  resulting  from  the  substitu¬ 
tions  for  83  and  82  in  terms  of  and  #1  as  in  [13].  Certain  qualitative  aspects  of  the 
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reduced  dynamics  can  still  be  explored,  sidestepping  analytic  difficulties.  For  instance  one 
can  investigate  relative  equilibria. 

Definition  6.1:  ze  £  M  is  a  relative  equilibrium  for  Xh  if 

*»(»(*.))  =  0-  (6-7) 

Remark:  Let  F*Xh  be  the  flow  of  Xh  on  M.  Then  2e  is  a  relative  equilibrium  iff  FjCjI(ze) 
is  a  stationary  motion,  i.e.  there  exists  £  £  Q  such  that 

fxhM  =  exp(tt)(ze). 

Theorem  6.2:  Let  J  be  an  Ad*  —  equivariant  momentum  mapping  on  M .  ze  £  M  is  a 
relative  equilibrium  of  Xh  iff  there  exists  a  £  £  Q  such  that  ze  is  a  critical  point  of 

ff(  =  H  ~  J(0  (6.8) 

where  J(£)  :  M  >  R  :  x  J(a:)(£). 

For  proof,  see  chapter  4  of  [l].  This  theorem  can  be  applied  to  a  simple  mechanical 
system  with  symmetry,  ( Q,K,V,G )  [19]. 

Theorem  6.3:(SmaZe)  For  simple  mechanical  system  with  symmetry  (Q,  K,  V,  G),  define 

V(  :  Q  R  :  q  ->  V(q)  -  \l<((Q(q),  ( Q(q ))  (6.9) 

for  each  £  £  Q .  Then  ze  =  (<jfe  ,  Pe)  G  T*  Q  is  a  relative  equilibrium  iff  qe  is  a  critical  point 
of  V$  for  some  £  £  Q  and  pe  =  Kb(£c}(qe)). 

Remark:  It  can  be  shown  that,  for  a  given  £  £  Q ,  has  the  symmetry, 

=  Vt(x)  (6.10) 

for  all  g  £  :=  {g  £  G\Adg£  =  £}.  If  G  =  S1  =  G  and  action  <3>  is  free  and  proper. 

Then  Q/G$  is  a  smooth  manifold  and  ^  :  Q  -»  Q/G$  is  a  submersion.  Thus  induces 
a  function  on  Q/G$  such  that 


O  7T  £ . 
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Recalling  that  V  =  0,  this  theorem  can  be  directly  applied  to  the  floating  four-bar 
linkage  system.  Rewriting  the  kinetic  energy  given  in  section  4  by  setting  the  origin  of 
inertial  space  at  the  system  of  center  of  mass,  we  have 

m  1 

T=-<u,3u>  (6.11) 

where  u>  =  (o>o ,  uj\  ) .  It  can  be  shown  that  the  infinitesimal  generator  on  Q  =  S1  x  S1  is 
Cq(q)  —  (1,1)T  [21].  Hence,  by  theorem  6.3,  (0*,/U)  is  a  relative  equilibrium  point  on 
T*Q  iff  9*  is  a  critical  point  of  function 

V<(«o,«>)  =  -(l,l)j(j) 

Applying  Proposition  4.2  and  the  definition  of  J  given  in  Theorem  4.1,  we  have 

^(^0,^1)  =  -eTJe  (6.12) 

where  e  =  (1  1  1  1)T .  From  Eq.(6.12),  we  observe  that  the  diagonal  terms,  which  include 
the  inertia  of  bars,  of  matrix  J  do  not  effect  the  positions  of  critical  points  of  function  . 
Since  the  elements  of  matrix  J  are  only  functions  of  0i—0o,  above  VJt  satisfies  (6.10)  for  all 
g  G  S1 .  It  follows  that  ^(#10)  =  ^(^Oj^i)-  Then,  the  critical  points  of  V^,  (0%,  9*),  will 
make  (9%  —  9(*)  to  be  the  critical  points  of  .  At  relative  equilibrium,  the  relative  angles 
between  bars  are  fixed  and  the  whole  system  rotates  around  the  system  center  of  mass 
with  constant  angular  velocity.  Unlike  the  planar  two-body  case,  the  relative  equilibrium 
shapes  depend  on  the  values  of  the  masses  and  the  lengths  of  bars.  From  the  expression 
of  matrices  J  and  the  relations  between  relative  angles  given  in  [13],  it  is  very  difficult  to 
find  critical  points  of  function  Rj  analytically.  However,  numerically  searching  for  critical 
points  is  easy  since  now  is  only  a  function  of  one  variable.  In  the  following,  we  give  an 
example  to  computing  the  relative  equilibria  by  applying  above  theorem. 

Example: 

An  assembly  in  Fig.  2  is  a  possible  structure  of  a  robot  arm.  The  parameters  are 
given  as  follows 

mo  =  3,  7ni  =  1,  m2  =  30,  m3  =  1; 
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Fig.  2 

d(0,3)  =  (-1.5,0),  d(0, 1)  =  (1.5, 0); 

d(l,0)  =  (-0.5,0),  d(l,2)  =  (0.5,0); 

d(2, 1)  =  (—6, 0),  d(2, 3)  =  (—3,0); 

d(3,2)  =  (-0.55,0),  d(3,0)  =  (0.55,0). 

Thus,  the  lengths  of  links  are 

Iq  =3,  W  —  I,  h  =  3,  h  —  1-1 

It  is  clear  that  s  +  l  <  p  +  q  and  s  =  /j  are  satisfied. 

The  graph  of  V, f  is  given  in  Fig.  3.  From  this  figure  we  can  see  two  critical  points 
appear  at  #10  =  7  deg  and  0\  0  =  186  deg .  The  shape  determined  by  these  two  angles  are 
shown  in  Fig.  4. 

19 


Fig.  3 


7.  CONCLUSIONS 

In  this  paper,  we  gave  sufficient  conditions  which  make  the  configuration  space  of 
a  closed  four-bar  mechanism  be  a  smooth  manifold  by  applying  the  classical  theory  of  a 
mechanisms  due  to  Grashof.  Under  one  of  them,  i.e.  s  +  l  <  p  +  q  we  derived  well  defined 
expression  of  kinetic  energy,  or  Lagrangian.  This  Lagrangian  is  invariant  under  the  action 
of  S 22(2) ,  the  rigid  motion  group  in  plane.  It  turns  out  that  the  four-bar  mechanism  is 
a  simple  mechanical  system  with  symmetry.  Applying  Poisson  reduction  we  obtained  the 
reduced  dynamics.  Furthermore,  by  using  Smale’s  theorem  on  relative  equilibria  we  found 
a  function  whose  critical  points  give  the  relative  equilibria  for  our  system.  An  example  of 
its  application  was  given. 
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